Frustrated Heisenberg Antiferromagnets between d = 2 and d = 3 
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The anitiferromagnetic Heisenberg model on a stacked triangular geometry with a finite number 
of layers is studied using Monte Carlo methods. A topological phase transition occurs at finite tem- 
perature for all film thicknesses. Our results indicate that topological excitations are important for 
a complete understanding of the critical properties of the model between two and three dimensions. 

PACS numbers: 75.10.Hk, 75.40.Cx, 75.40.Mg 
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The behavior of the Heisenberg antiferromagnet on tri- 
angular geometries is quite unusual. The ground state or- 
der parameter has SO (3) symmetry which allows for a Z2 
topological defect. In d = 2, there is a finite temperature 
vortex unbinding transition which is purely topological in 
characteri2i2i2i£iii. Both above and below this tempera- 
ture the spin-spin correlation function decays exponen- 
tially with distance. This behavior is quite different from 
the Kosterlitz-Thouless transition 7,8 which occurs in two- 
dimensional XY models. In the latter case the correla- 
tions are exponential at high T and decay with a power 
law at all temperatures below the transition. There is no 
long range order but there is a finite spin wave stiffness 
in the low temperature phase. In the Heisenberg anti- 
ferromagnet, both the sublattice magnetizations and the 
spin wave stiffness are zero at all finite temperatures. On 
the other hand the vorticity modulus is nonzero at low 
T and drops to zero at the transition. Hence the phase 
transition is similar to the Kosterlitz-Thouless transition 
in that vortices are involved but it is different in that 
the correlations decay exponentially at all finite temper- 
atures. 

Field theoretic studies of this model in terms of the 
nonlinear sigma model (NLa) predict that there is only 
a transition at T c = in two dimensions. However a 
fixed-dimension perturbative approachSiifi in d = 2 us- 
ing the corresponding Landau-Ginzburg- Wilson (LGW) 
Hamiltonian has located a fixed point which may de- 
scribe a topological phase transition. In d — 3 both the- 
ory and experiment indicate that there is a conventional 
phase transition with the appearance of an order param- 
eter below the critical temperature. However there is a 
great deal of debate about the nature of this transition. 
Monte Carlo studies and perturbative field theory calcu- 
lations in fixed dimension 11 indicate a continuous tran- 
sition which belongs to a new chiral universality class 
whereas non-perturbative RG (NPRG) studies indicate 
a phase transition which is very weakly first order with 
effective exponents^^*^. A very recent numerical study 
of the RG flow in the LGW Hamiltonian by ItakuraM 
also suggests a possible weak first order transition for 
the Heisenberg spin model but with much stronger ev- 
idence for a first order transition in the case of XY 



spinsi^i^. The experimental studies of these three di- 
mensional stacked triangular systems find that they ex- 
hibit a continuous phase transition with a set of critical 
exponents that do not belong to any of known universal- 
ity classes . 

In this work we explore the crossover of the d = 2 
behavior to d = 3. We study the nature of the ordering 
process of a layered triangular system having dimensions 
L x L x H with the number of stacked triangular layers 
ranging from H = 2 to H — 24 and the linear size of each 
layer L ranging from 18 to 120 by means of Monte Carlo 
simulations. 

We consider an isotropic Heisenberg model of classical 
spins interacting via nearest neighbour exchange on a lay- 
ered triangular lattice described by the following Hamil- 
tonian 



H — JijSi ■ Sj 



(1) 



Within each triangular layer the spins interact antifer- 
omagnetically Jy < 0. Between adjacent layers the in- 
teraction can be taken to be positive or negative with 
no essential difference since there is no frustration asso- 
ciated with these interactions and hence there is a sym- 
metry with respect to the sign. We shall choose it to be 
ferromagnetic Jy > and set | Jy | = J = 1 for all inter- 
actions. Periodic boundary conditions are applied in the 
plane of each triangular layer where, in order to preserve 
the 3-sublattice structure, it is important to choose L to 
be a multiple of 3. The finite layer thickness, H, can have 
any integer value since we have chosen the inter-planar 
interactions to be ferromagnetic. The top and bottom 
surface layers are subject to free boundary conditions. 

The spin stiffnessi*^, or helicity modulus, is a mea- 
sure of the increase in free energy associated with a twist- 
ing of the order parameter in spin space by imposing a 
gradient of the twist angle about some axis n in spin 
space along some direction u in the lattice. Diagonal 
elements of the spin wave stiffness tensor can be calcu- 
lated by choosing an orthogonal triad for the directions 
of three principal axis in spin space. The symmetry of 
the ground state suggests that two of the principal axes 
correspond to two perpendicular directions and J_2 
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in the spin plane and that the third is perpendicular to 
this plane This third axis is conveniently chosen to point 
along the average chirality direction K. 

Another quantity of interest when topological proper- 
ties of the system are to be considered is the vorticityi 
The vorticity is a measure of the response of the spin 
system to an imposed twist about an given axis h in spin 
space along a closed path that encloses a vortex core. 
The vorticity 14 contains a contribution due to the vor- 
tex core as well as part which is proportional to ln(L/a). 
This can be written as 

Vn = C A + v il \n(L/a). (2) 

where Cn is a temperature dependent constant describing 
the core and Vn is called the the vorticity modulus. It 
plays a similar role to the spin stiffness and vanishes at 
a phase transition. 

Azaria et alm^ have made detailed predictions for the 
dependence of the stiffnesses on the linear size of the sys- 
tem for single layers based on the continuum limit of this 
model. Their results for a nonlinear sigma model (NLct) 
using renormalization group (RG) techniques showed 
that the spin wave stiffness of the triangular antiferro- 
magnet is a nontrivial function of In L at low tempera- 
tures where L is the linear dimension of the triangular 
lattice. At any finite temperature T, the stiffnesses are 
zero on length scales large compared with the correla- 
tion length. However, on length scales 1 <C L <C £, the 
stiffnesses are nonvanishing at low T. 

In order to make a comparison of our calculated stiff- 
nesses with the predictions of the RG equations, we first 
calculate the stiffnesses for linear size L = 18 so that 
the condition H <C L is satisfied. Using the values at 
this length scale, we numerically integrate the expres- 
sions of Azaria et. al. with these initial values to predict 
the behavior at larger length scales. We then compare 
the predicted curves with the values obtained directly 
using Monte Carlo methods at the relatively low tem- 
perature T = 0.2 for different layer thicknesses. The de- 
pendence of the average of the three principal stiffnesses, 
p = 5 En=i Pn> 0n m ^ * s sri0wn m figure n for various 
values of H . In this figure the solid lines are the low 
temperature T = 0.2 RG predictions and the points are 
the Monte Carlo results. The data points agree very well 
with the predictions of the NLc model at this low tem- 
perature. The stiffnesses decrease proportional to InL 
and the results indicate that there is no stiffness at large 
length scales. Hence it would appear that the finite layer 
systems behave in the same way as the single layer system 
with no indication of a finite temperature phase transi- 
tion in the stiffnesses. However, as the the number of lay- 
ers increases the slopes of the lines in figure ^ decrease 
indicating that the stiffnesses should approach a finite 
value at large length scales as the system becomes fully 
three dimensional. 

The excellent agreement of the Monte Carlo results 
and the RG predictions above indicate that topological 
excitations are not important in this system at low T. 
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FIG. 1: Average spin stiffness plotted as a function of InL for 
different layer thicknesses H at T — 0.2. The solid lines are 
the low temperature RG predictions. 
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FIG. 2: Vorticity data for H = 2 and different linear sizes 
L — 18, 30, 42, 60, 72 as function of temperature. 

However, at higher temperatures they may play an im- 
portant role. These topological degrees of freedom are 
not included in the RG analysis. In order to study the 
vortex degrees of freedom directly, we have also calcu- 
lated the response of the system to the presence of a sin- 
gle virtual vortex. We calculate this vorticity response 
for a fixed number of layers H and for various linear 
sizes L. In figure |21 the temperature dependence of the 
vorticity data for H = 2 is shown. The curves for dif- 
ferent sizes L all cross at a common temperature point 
near T = 0.51. This behavior suggests that there is a free 
energy cost to creating isolated vortices in the system be- 
low this temperature but that they should spontaneously 
appear above it. At any value of T we can study the size 
dependence of the vorticity using equation @ to extract 
the coefficient of the In L term which we call a vorticity 
modulus. 

The vorticity moduli have been determined for systems 
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FIG. 3: Temperature dependence of H times the average vor- 
ticity moduli for H = 2, 3, 4, 6, 8 and 20. The solid line rep- 
resents the ratio Hv/T = —. 
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FIG. 4: T C (H) dependence on 1/H is indicated by the solid 
circles and the temperatures of the specific heat maxima are 
indicated by open circles. The dotted line represents a fit to 
the prediction of FSS theory given in and the solid line is 
the same fit with H± = 0. 



having a different number of layers H. In each case, the 
vorticity moduli abruptly drop to zero at the same tem- 
perature where the vorticity curves cross. Figure [3] shows 
the average vorticity moduli v = | J2n=i v n multiplied 
by the layer thickness H as a function of temperature for 
several layer thicknesses. The solid line represents the ra- 
tio Hv/T = — which intersects the corresponding layer 
vorticity modulus at the temperature where the raw vor- 
ticity curves cross. This suggests that the average vortic- 
ity modulus exhibits a universal jump , v/T c = 1/(ttH), 
at the transition which decreases as the number of layers 
H increases. 

Figure|3]shows a plot of T c as well as the temperature of 
the specific heat maximum versus the inverse number of 
layers . As -tt — > the system becomes three dimensional 
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FIG. 5: Average vorticity v(z) for H = 14 for each individual 
layer at T = 0.8. The dotted line is the average over all layers. 

with the critical temperature^ equal to T c (3d) = 0.958. 
The specific heat maxima lie above the temperatures 
where the vorticity modulii drop to zero but the two tem- 
peratures merge together as 1/H — ► 0. 

The data points can be fit very well by a slightly mod- 
ified form of the usual finite scaling (FSS) 

T C (H) = T c (3d) (l C \ — ) (3) 
V (H + Hi) « J 

where both C\ and H\ are constants. The correlation 
length exponent v in the above expression is the value 
for the three dimensional cassia and equals v= 0.59. At- 
tempted fits of the data in figure 0] (solid line) using this 
form with H \ = are are only possible for large values 
of H > 8 and there are strong deviations for smaller val- 
ues of H. The dotted line in figure 0] corresponds to the 
values C\ — 7.28 and H\ — 3.06. This modified form 
provides an excellent fit over the entire range of layer 
thicknesses. 

The quantity H+Hi behaves as an effective layer width 
which may be related to the boundary conditions at the 
upper and lower surfaces. Figure show the average vor- 
ticity modulus v(z) for each individual layer for the case 
of H = 14 . The dotted horizontal curve is the average 
over all layers. The deviations of v(z) from the average 
are larger for thin films than for thicker films. For small 
values of H the effect of surface layers is significant and 
in the vicinity of the critical point the thickness H is not 
the only relevant length scale. It is being modified by the 
additional constant H\ which is related to how quickly 
v{z) changes near the surface due to the free boundary 
conditions. 

Modifying finite scaling theory to include the effective 
width for small H < 8 predicts that the vorticity moduli 
v should scale with respect to the layer thickness H as 
follows 

{H + H 1 )v = f{{H + H 1 ) l ^\t\). (4) 
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FIG. 6: FSS scaling plot of the average vorticity moduli v. 



where v = 0.59 is the three dimensional correla- 
tion length exponent and t is the reduced temperature 
(T c (3d) - T C (H))/T c (3d). FigureEJshows a plot of H eff v 

in terms of the variable H^L v At\. The vorticity moduli 
curves for different thicknesses collapse onto a universal 
curve. 

The results above are similar to those found in studies 
of 4 He superfluid nmsi2i22i2I. In these systems the or- 
der parameter is a complex number which has the same 
symmetry as the ferromagnetic XY model. The super- 
fluid transition in d = 2 can be described in terms of 
the Kosterlitz-Thouless transition where vortices destroy 
the ordered phase. Ambegaokar et. alJ& argued that the 
KT theory can be extended to films of finite thickness and 
that the transition should continue to have the 2d charac- 
ter provided the thickness is not too large. Schultka and 
Manousakis— have studied the same problem numeri- 
cally and found that the results could indeed be described 
by the KT theory using finite scaling theory. However, 
a detailed fit for a small number of layers required the 



introduction of an effective width for the filmsiSi Our re- 
sults for the frustrated Hciscnberg model are quite sim- 
ilar but there are some important differences. In the 
XY model, the transition is accompanied by a jump in 
the stiffness with a power law decay of the correlations 
below this temperature. This behavior perists for films 
of arbitrary thickness until the 3c? bulk limit is reached 
where the transition becomes the usual A transition. In 
our case, the d = 2 behavior observed^^^ for a single 
layer also persists as the number of layers increases. The 
stiffness is zero at all finite T indicating that the spin cor- 
relations decay exponentially at all temperatures. How- 
ever, the vorticity modulus exhibits a universal jump at 
a finite temperature indicating that a topological phase 
transition occurs involving vortices. The size of the jump 
varies as l/H and hence approaches zero in the limit 
of d = 3. In addition, the stiffness also approaches a 
nonzero value at temperatures below the critical tem- 
perature in this same limit. This indicates that the 2d 
behavior approaches the 3d behavior in a continuous way 
and that topological excitations play an important role 
between the two limits. The actual role that they play 
in determining the nature of the transition in 3d is still 
not clear. It is possible that the Heisenberg model has 
a weak first order transition as predicted by the NPRG 
field theories and that the critical exponents are only ef- 
fective exponents. However, the evidence for first order 
behavior is much clearer for the XY model and the dif- 
ferences between Heisenberg and XY behavior may in 
fact be due to the different types of topological defects. 
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